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1. Introduction 

Strongly correlated electronic systems are currently intensely studied in relation with high T c 
superconductivity. Recently there has been a renewed interest in the one-dimensional t-J model as 
an integrable low dimensional version of a strongly correlated electronic system. The t-J model was 
proposed by F.C. Zhang and T.M. RiceW. P.W. Anderson claimed that two-dimensional systems 
may share features of one-dimensional systems' 2 ', which could imply a certain relevance of some 
results obtained for the one-dimensional t-J model to high T c superconductivity. 

The model describes electrons on a one-dimensional lattice with a Hamiltonian that includes 
nearest neighbour hopping (t) and nearest neighbour spin exchange and charge interactions (J). 
The Hilbert space of the model is constrained to exclude double occupancy of any single site, which 
corresponds to an infinite on-site repulsion. 

Electrons on a lattice are described by operators Cj a , j = 1, . . . , L, a = ±1, where L is the total 
number of lattice sites. These are canonical Fermi operators with anti-commutation relations given 
by {c\ a ,Cj T } = S i jS CF T . The state |0) (the Fock vacuum) satisfies c j cr |0) = 0. Due to the constraint 
excluding double occupancy there are three possible electronic states at a given lattice site i 

io) , 1^=410), 1 1), =4-ii°) • a- 1 ) 

By n i a — c\ a c i tJ we denote the number operator for electrons with spin a on site i and we write 
»i - + The spin-operators S = £^ =1 Sj, = £^ =1 Sj, S* = £^ =1 S?, 

Sj = Cj t iCj t -i , Sj — c]^iCj i , Sj — — — n j,-i) (1-2) 

form an su(2) algebra and they commute with the hamiltonian that we consider below. [We shall 
always give local expressions Oj for symmetry generators, implying that the global ones are obtained 
as = Y^j=i The Hamiltonian on a lattice of L sites is given by the following expression : 

L 

3 = 1 a=±l 

3 = 1 V 

The projectors (1 — rij _ a ) in the kinetic energy terms ensure that H acts within the constrained 
Hilbert space, i.e. does not create states with double occupancy. A equivalent expression for H is 

H = E [~ tV E ( C LS + 1, CT + h.c) V + J (Sj ■ S j+1 \nj n j+1 ) | , (1.4) 

where V is the projector on the subspace of non-doubly occupied states. It was shown in [3] that 
the model described by (1.3) is integrable and can be solved by coordinate space Bethe Ansatz. In 
addition it is possible to map it onto the integrable quantum lattice gas of Lai and Sutherland I 4 ' 5 ! . 

The thermodynamics of the model were treated in [3,6] and the ground state properties and 
excitation spectrum were investigated in [7-9]. 
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As the number operator for electrons N = S CT =±i n j,a commutes with H, we can add 

a term 2 N — L to the hamiltonian without changing the set of eigenvectors. Physically this just 
amounts to a shift of the chemical potential. For the special value 2 J = 2t = 2 the shifted hamiltonian 
now exhibits a number of interesting properties: 

Firstly it is supersymmetric, i.e. it commutes with all nine generators of the superalgebra w(l|2)[ 6,7 l. 
Secondly it can be written as a graded permutation operator! 6 ' 7 ]: 

L 

H susy = H + 2N-L = -Y / ^ +1 

(1.5) 

J - 1 susy 

i=i 

The operator permutes the three possible configurations (1.1) between the sites j and j + 1, 

picking up a minus sign if both of the permuted configurations are fermionic, i.e. 

n^+^o),. x |o) i+1 = io),. x |o) i+1 

n« +1 10)^ x = 1^. x |0> j+1 (1.6) 

U^+ 1 \r) j x \a) j+1 = -\a) j x |r) j+1 , a,r =T,| . 

It is clear that this form of interaction conserves the individual numbers and of electrons 
with spin up and spin down, and due to the constraint of no double occupancy the number N h of 
empty sites (or "holes") is also conserved. 

The outline of the paper is as follows : 
In section 2 we give a discussion of the supersymmetry of the model and express the hamiltonian 
(1.5) in a way most suitable for the analysis of section 5. In section 3 we perform a detailed con- 
struction of the Algebraic Bethe Ansatz of the model. We derive three different forms for the Bcthc 
Ansatz Equations (BAE) and the eigenvalues of the transfer matrix. Two of these solutions have 
previously been obtained by various authors' 3-5,10 ', whereas our third solution, presented in section 
4, is new. Our expression for the BAE seems to be particularly simple and we hope that it will be 
useful in clarifying the physical features of the model. The direct physical consequences of our new 
solution (like the structure of the ground states and classification of excitations and correlations) are 
currently under investigation' 11 ! and will presented in a future publication. The graded Quantum 
Inverse Scattering Method (QISM), discussed in section 3, enables us to obtain expressions for (an 
infinite number of) higher conservation laws at the quantum level. These conserved charges are 
of interest, because physical interactions, although of short range, are not generally well approxi- 
mated by interactions involving only nearest neighbours. The charges under consideration involve 
interactions of longer range (next nearest neighbours, next next nearest etc.) and can be added to 
the hamiltonian without spoiling the integrability of the model. Thus it is possible to construct 
integrable models with longer range interactions by using higher conservation laws' 12 '. The first 

2 The case —J = 2t can be obtained from J = 2t via the transformation 

Cj - (-1)^ , c] -» i-iy c ] 
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nontrivial higher integral of motion is for example given by the expression 

L 

*;=i 

where H k ' k+1 is the density of the hamiltonian defined in (1.5). Section 5 is devoted to the derivation 
of explicit formulas for higher conservation laws. 



2. Supersymmetry of the t-J model 



For J — 2t the t-J model exhibits a (global) u(l\2) invariance on the constrained Hilbert space. 
In the literature this symmetry algebra has frequently been denoted by spl(2,l)W. The relation 
between spl(2, 1) and sw(l|2) (neglecting the trivial u(l) factor for the time being) is as follows : 
The algebra sw(l|2) is a real form of the complex Lie-superalgebra sZ(l,2; C), whereas spl(2, 1) is an 
equivalent notation for sl(l, 2; C). For our purposes it is more convenient to work with a real base 
field, so that we will work in a representation of w(l|2). The generators of the u(l|2) algebra are 
given by I 13 l (we write O = X^t=i ®j) 



r - <?t 



C ),l C 3 -1 










J j;6 






1, 


-i) 






J 3,8 


= ( 1 -« J -,-i) 








= ( 1 -»i,-i) 







(2.1) 



The operators 5, S\ S z , Q 1 , Q\, Q_i, Q^i, T generate the su(l|2) subalgebra of w(l|2). 
In the fundamental representation there exists an invariant, nondegenerate bilinear form K af3 
on u(l|2), which is given as the supertrace over two generators I 14 l 

K aP ={K^y 1 = S br{J j<a J jtP ) (2.2) 



and which is explicitly computed in Appendix A. Note that the nondegeneracy of K a/3 is a feature 
of the fundamental representation and docs not generally hold for other representations, because 
u(l|2) is not scmisimple. For later use we define the structure constants of w(l|2) as 

I4ai4^} : = Jk.aJk.f) ~ (~l) e " e,3 <4 :/ 3<4 : a = /a/3 7 -4, 7 > ( 2 '3) 

where e a are the Grassmann parities of the generators J k a (i.e. e = 1 for the fermionic generators 
J 4 . . . J 7 and e = for the rest). 

The Hamiltonian of the t-J model on the constrained Hilbert space can now be expressed in 
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terms of the densities J ka as 

L L 

H — FPJ+ 1 — — K af3 T T 
11 susy — / j 11 susy ~ / , " J j,a j+1,/3 

3=1 3=1 
L 

= -E E Qj,*Q>+i,* + Qj+i,*Qj,* ( 2 - 4 ) 

j = l <T=±1 

L 

- ^2 -2SjS* +1 - S^S j+1 - SjS* j+1 + 2TjT j+1 - 

3 = 1 

where we have used (2.1) and the explicit expression for K af3 given in Appendix A. In this form the 
global u(l\2) invariance of the hamiltonian is easily confirmed 

[H,J a ]=0 a=l,...9 . (2.5) 

The "manifestly supersymmetric" expression for the hamiltonian (2.4) in terms of the Killing form 
will enable us to express the higher conservation laws we derive in section 4, in an u(l|2) invariant 
way. 

3. Graded Quantum Inverse Scattering Method 

In this section we construct the enveloping spin model of the hamiltonian of the one-dimensional 
supersymmetric t-J model, using the Quantum Inverse Scattering Method (QISM). Due to the fact 
that we are dealing with a supersymmetric theory it is necessary to modify the QISM along the lines 
discussed in [10,15]. Below we give a summary of the "graded" version of the QISM. We start with 
an i?-matrix, obeying a graded Yang Baxter equation, and construct from it a "fundamental" spin 
model (i.e. the L-operator is constructed directly from the i?-matrix). We then show that the trace 
identities of the corresponding transfer matrix give rise to the hamiltonian of the t-J model. Finally 
we construct a set of simultaneous eigenstates of the transfer matrix and the hamiltonian, using a 
nested Algebraic Bethe Ansatz 3 (NABA)t 10 ' 19 ~ 21 l . Due to the grading there exist three choices for the 
i?-matrix, all of them describing the same physical system, but leading to different (yet equivalent as 
shown in Appendix C) forms of the NABA. Two of these possibilities of performing a Bethe Ansatz 
analysis are equivalent to the coordinate space Bethe Ansatz solutions of Lai and Sutherland. We 
reproduce their respective periodic boundary conditions and expressions for energy eigenvalues of 
the hamiltonian. 

3.1. Yang Baxter Equation 

Consider the graded linear space y( m l") = V m V n , where m and n denote the dimensions of 
the "even" (V m ) and "odd" (V n ) parts, and © denotes the direct sum. Let {e 1 , . . . , e m+n } be a 
basis of v^ m+n \ such that {e 1; . . . , e m } is a basis of V m and {e m+1 , . . . , e m+n } is a basis of V n . The 

3 For a general introduction to the Algebraic Bethe Ansatz we refer to the publications (16-18). 
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Grassmann parities of the basis vectors are given by e x = . . . = e rn = and e rn+1 = . . . = e m+n = 1. 

Linear operators on T/( m l") can be represented in block form (M € End(V^ m ^)) 

(A B\ (A 0\ /0 B\ 

M = ( , e = , e = 1 . (3.1) 

\C Dj \0 D J \C J 

The supertrace is defined as 

str(M) = tr(A) - tr(D) , (3.2) 

where the traces on the r.h.s. are the usual operator traces in V m and V n . We now define the graded 
tensor product space T/( m l™) ® y( m \ n ) j n terms of its basis vectors {e a ® e b \a,b = 1, . . . ,m + n} as 
follows 

v®w = (e a v a )®{e b w b ) = {e a ®e b )v a w h {-iy^ f - b . (3.3) 

Compared to the "ordinary" tensor product the additional factor (— l) 6 "" 66 occurs, which originates 

in passing v a past e b . The action of the right linear operator F®G on the vector v ® w in T/( m l") 
y(m|n) is given by 

(F ® G) {v ® w) = F{v) ® G{w) . (3.4) 
Therefore its matrix elements are of the form 

{F®G)f d = F ab G cd {-lf^ + ^ ■ (3-5) 
The identity operator in \/( m l") ® yi m \ n ) j s given by I^i, 1 = ^a n 6 n ^o,6, an( i t nc matrix II that 



permutes the individual linear spaces in the tensor product space, is of the form 

Yl(v®w) = (w®v), (U)ll b b l=S aib2 5 a2bi (-iy^ . (3.6) 

The physical relevance of the above construction is as follows : If we consider a lattice gas of m 
species of bosons and n species of fermions, then \/( m l") denotes the space of configurations at 
every site of the lattice. For the example of the t-J model we have m = 1, n — 2 and the three 
allowed configurations are given by (1.1). The tensor product space l/W 2 ) ®V^ 2 ^ describes two 
neighbouring sites, and owing to the fermionic nature of some of the configurations, the tensor 
product has to carry a grading. II permutes configurations on neighbouring sites, and we pick up a 
minus sign if we permute two fermions. 

A matrix R(X) (depending on a spectral parameter A) is said to fulfill a graded Yang-Baxter- 
equation, if the following identity on V {m ^ ® V (m ^ ® V^ m ^ holds 

{I ® R{\ - H)) {R(X) ® I) {I ® R(n)) = {R(n) ® I) (I ® R(X)) {R(X - H) ® I) . (3.7) 

In components this identity reads 

R(X riZZ R(X) a c f d2 R(n)%% = R^ZZ RW C a 2 X J2(A - . (3.8) 

Note that despite the fact that the tensor product in (3.7) carries a grading, there are no extra signs 
in (3.8) compared to the nongraded case. It is easily checked that the R-matrix 

R(X)=b(X)I + a(X)U 

A i (3-9) 

fulfills equation (3.8). 
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3.2. Construction of the Transfer matrix 



By multiplying (3.8) by 11^ lllfal from the left one can derive the equation 

R(\-^lZ (HR(\)) C /X (ni?( M ))^ (-i)^ 2 (%+%) = 

(3.10) 

(ni2(A0)j£ (nfl(A))^ J?(A-^fc (-l)^^ . 

In matrix notation (3.10) reads 

i? 12 (A - /i) ([n 13 i? 13 (A)] ® [u 23 rm]) = ([n 13 i? 13 ( M )] ® [n 23 i? 23 (A)]) r 12 (\ - n) , (3.11) 

where the indices 1, 2, 3 indicate in which of the spaces V^-" 1 ^ in the tensor product space y( m l") <gi 
y(m\n) ^ y{m\n) ^ e matrices act nontrivially. The tensor product in (3.11) is between the spaces 
1 and 2. We now call the third space "quantum space" and the first two spaces "matrix spaces". 
The physical interpretation of the quantum space is as the Hilbert space over a single site of a 
one-dimensional lattice. We now consider the situation, where intertwining relations of the type 
(3.11) hold for all sites of a lattice of length L. The quantum space index "3" now gets replaced by 
an index labelling the number of the site. We define the L-operator (on site n) as a linear operator 
on TL n ® V^iatrix (where TL n ~ V^ m ^ is the Hilbert space over the n th site, and V^}^\ x is a matrix 
space) 

- n~i2(A)$ = (fe(A)H + a{X)I)% . (3.12) 

L n is a quantum operator valued (m + n) x (to + n) matrix, with quantum operators acting non- 
trivially in the n th quantum space (of the direct product Hilbert space over the complete lattice 
®j =1 Ti.j). The greek indices are the "quantum indices" and the roman indices are the "matrix 
indices". Equation (3.11) for the n th quantum space can now be rewritten as the operator equation 

R{X-n)Zc\ L n {X%X LMf* (-1)MS+S) = 

(3.13) 

In matrix notation (3.13) takes the form 

R(X - ft) (L n (X) ® L n (u)) = {L n {lA ® L n {X)) R(X - M ) . (3.14) 

Here the graded tensor product is between the two matrix spaces and R only acts in the matrix 
spaces. The intertwining relation (3.14) enables us to construct an integrable spin model as follows 

We first define the monodromy matrix T L {X) as the matrix product over the L-operators on all sites 
of the lattice, i.e. 

T L (X)=L L (X)L L _ 1 (X)...L 1 (X) 
((T L (X)r»)c . .. ^JA)^ [ L i _ 1 (A ) :';-_ i ...L 1 (A)5 i x (3 . 15) 

Pi—Pl 
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T L (A) is a quantum operator valued (m + n) x (m + n) matrix that acts nontrivially in the graded 
tensor product of all quantum spaces of the lattice and by construction fullfills the same intertwining 
relation as the L-operators (as can be proven by induction over the length of the lattice) 

R(X - M ) (T L (A) ® T L (p)) = (T L (p) ® T L {\)) R(X - fi) . (3.16) 

The transfer matrix r(A) of the integrable spin model is now given as the matrix supertrace of the 
mondromy matrix 

m+n 

r(A) = str(T L (X)) = £ (-l)^(T L (A)) aa . (3.17) 

a=l 

As a consequence of (3.16) transfer matrices with different spectral parameters commute. This 
condition implies that the transfer matrix is the generating functional of the hamiltonian and of an 
infinite number of "higher" conservation laws of the model. 



3.3. Trace Identities 

Taking logarithmic derivatives of the transfer matrix at a special value of the spectral parameter, 
one can generate higher conservation laws! 16 !. For our specific case at hand, i.e. the i?-matrix (3.9), 
the corresponding hamiltonian is obtained by taking the first logarithmic derivative at zero spectral 
parameter 

L 

= -V(n fe < fe+1 - 1). (3.18) 



dln(r(X)) 



*=0 fe=l 



The proof of this identity can be carried out in the same way as for the ungraded case, the main 
difference being the grading of the tensor product of the quantum spaces (see (3.15)). By shifting 
the energy eigenvalues by a constant we obtain the expression (1.5) for the hamiltonian of the t-J 
model 



dln(r(X)) 

U susy 1 



- L = H {2) - L , (3.19) 

A=0 



if we choose our underlying graded vector space to have signature (1,2), i.e. to have a basis with 
two fermionic and one bosonic state. This shows that the transfer matrix constructed from the 
L-operator (3.12) and i?-matrix (3.9) is indeed the correct transfer matrix for the one-dimensional 
supersymmetric t-J model. Higher conservation laws are obtained as the coefficients of the power 
series 

OO . 

ln(r(X)(r(0))-')=Y / ^H (k+1) . (3.20) 

k=l 

There exists however a simpler method for the construction of higher integrals of motion than taking 
logarithmic derivatives, which we will discuss in section 5. 
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3.4. Algebraic Bethe Ansatz with a bosonic background (FFB grading) 
Lai solution 



Due to the constraint of no double occupancy there are three different configurations per site 
for the t-J model. Thus the Hilbert space at the k th site of the lattice is isomorphic to C 3 and 
is spanned by the three vectors e 1 = ( 1 0) T , e 2 = (0 1 0) T and e 3 — (0 1 ) T . In this 
section we consider a grading such that e 1 and e 2 are fermionic (representing spin down/spin up 
electrons respectively) and e 3 is bosonic (empty site) . In terms of the Grassmann parities this means 
that c-l = e 2 = 1 and e 3 = 0. We pick the reference state in the k th quantum space |0) fe , and the 
vacuum |0) of the complete lattice of L sites, to be purely bosonic, i.e. 



|0>„ = 



|0) = ®Li|0> r 



(3.21) 



This choice of grading implies that R(fi) = b(fi)I + a(fi)H is given by the following expression 



/6(A)-a(A) 




























6(A) 





-a(A) 























6(A) 











a(A) 











-a(A) 





6(A) 





























6(A)-o(A) 





























6(A) 





a(A) 











a(A) 











6(A) 























a(A) 





6(A) 





























1 



R(X) = 



\ 

The .L-operator is defined by (3.12) and is of the form 

/a(X)-b(X)e^ -b{X)ef 



(3.22) 



V 



-b(\y n 2 

b(\)el? 



b(X)ef 
b(X)ef 



\ 



(3.23) 



a(X)-b(X)ef 

b(X)el 3 a(X) + b(X)el s J 

where e° b are quantum operators in the n th quantum space with matrix representation (e° b ) Q/3 = 
5 aa 8 b p. The monodromy matrix (3.15) is a quantum operator valued 3x3 matrix, which we represent 



as 



T L (X)=L L (X)L L _ 1 (X)...L 1 (X) 
The transfer matrix is then given as 



/A u {\) A 12 (X) B^xy 
A 2l (X) A 22 (X) B 2 {X) 
\ C\(X) C 2 (X) D(X) 



(3.24) 



- str{T L (n)) = —A n (fi) - A 22 (fi) + D(n) 

The action of L k (X) on the reference state on the k th site is 

/ a(X) 

L k (X)\0) k = 



(3.25) 







\&(A)4 3 b(X)ef 



a(X) | |0) fc 



(3.26) 
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Using (3.15) and (3.26) we determine the action of the monodromy matrix on the reference state to 
be 

/(a(X)) L 0\ 
T L (X)\0) = (a(X)) L |0) . (3.27) 
V C 1 (\) C 2 (A) I J 

We will now construct a set of eigenstates of the transfer matrix using the technique of the NABA. 
Inspection of (3.27) reveals that C 1 (X) and C 2 (A) are creation operators (of odd Grassmann parity) 
with respect to our choice of reference state. This observation leads us to the following Ansatz for 
the eigenstates of t(/j) 

|A 17 . . . , X n \F) = C Bi (A x ) C a2 (A 2 ) . . . C 0b (AJ |0) , ( 3 .28) 

where the indices a,j run over the values 1,2, and F a ^--- a i is a function of the spectral parameters 
Xy The action of the transfer matrix on states of the form (3.28) is determined by (3.27) and the 
intertwining relations (3.16). The components of the intertwining relations (3.16) needed for the 
construction of the NABA are 

AM C C (X) = (-l)Vp r ^~_ A ^ g C p (X) A a M + fc^|c b ( M )A ac (A) 

DQ>) CJX) = _L_ C c( A) DO.) - ^CMD(X) ^ 
C 0l (A 1 )C 02 (A 2 ) = r(A 1 -A 2 )^ C &2 (A 2 ) C^) , 

where 

K^)cd = K^ab^cd - a{n)S ad 8 bc 

Here II^ 1 )"^ = — <5 ad S bc is the 4x4 permutation matrix corresponding to the grading e 1 = e 2 = 1. 
r(/j) can be seen to fulfill a (graded) Yang-Baxter equation on its own 

r(A - »)Z% r{X) a c ± r(») d c f h l = r( M )2£ r ( A )a't ~ ^1 > ( 3 ' 31 ) 



(3.30) 



and can be identified with the i?-matrix of a fundamental spin model describing two species of 
fcrmions. Using (3.29) we find that the diagonal elements of the monodromy matrix act on the 
states (3.28) as follows 



n , 

D( M )|A 1; . . . , X n \F) = ft — -|A X , . . . , X n \F) 



n b b n 

+EW 1 "' n c^Hc^mf^ 

f— ' V / a 1 ...a n " , 3 

k=l n 7 = 1 



(3.32) 



3- 



(A 11 ( M ) + A 22 ( M ))|A 1 ,...,A„|F) 



n 1 n 

n n 

fe=l j=l 



(3.33) 



where 

r (l) (M) V;^ =str(T (l) (M)) 

= str(L^ (fi - XjL^Ui, Vi) • • • 4 15 (M - A 2 )i 1 1) (a* - AJ) , 

and 

(3.34) 

4 1) (A) = 6(A)nW+a(A)/( 1 ) 

/a(A)-6(A) e i 1 -6(A) eg 1 \ 

= n«r(A) = 

V -6(A) e^ 2 a(A)-6(A) e f/ 

L^ 1 ) and r(/ti) can be interpreted as L-operator and i?-matrix of a fundamental spin model (r fullfills 
the Yang-Baxter equation (3.31)), describing two species of fermions. Hence Tn\fj,) and r^(fi) are 
the monodromy matrix and transfer matrix of the corresponding inhomogeneous model. Inspection 
of (3.32) and (3.33) together with (3.25) shows that the eigenvalue condition 

T(n)\\ 1 ,...,\ n \F) = v(f,,{\ j },F) |A 1 ,...,A n |F> (3.35) 

leads to the requirements that F ought to be an eigenvector of the "nested" transfer matrix t^^(/x), 
and that the "unwanted terms" cancel, i.e. 

(-(AfeAvi; + (h)a\-X) F a - ai = • (3-36) 

The relative sign in (3.36) is due to the supertrace in (3.25) and (3.35). The quantities A fc and 
A k are computed in Appendix A. Using their explicit expressions in (3.36) we obtain the following 
conditions on the spectral parameters A ■ and coefficients F, which are necessary for (3.35) to hold 

K A fc))- L f[ ^-"a'! Fb - bl = > k = l,...,n . (3.37) 

Ijik 

This completes the first step of the NABA. In the next step we will now solve the nesting. The 
condition that F ought to be an eigenvector of T^(fi) requires the diagonalisation ofr^(fi), which 
can be carried out by a second, "nested" Bethc Ansatz. From (3.31) and (3.34) the following 
intertwining relation is easily derived 

r(A - M ) (t«(A) ®tW(»)) = (tW^ ® TW(A)) r(A - . (3.38) 

If we write 



(3.39) 



then (3.38) and (3.30) imply that 

ft(A — (Jbj dyfj, — A J 

CW(A)C«(m) = C (1) (a^ (1) (A) • 
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As the reference state for the nesting we pick 

l )l 1) = (J) . |0>«=^ = i|0)i 1) • (3.41) 

The action of the nested monodromy matrix Tn\n)on the reference state 0}^ is determined by 
(3.34) and we find 



AW( A1 )|0)W = n«(M-A,)|0)( 1 ) 



DWM\0)W = fl (a(A* - A,) - b(, A,)) |0>« = f[ ^-^y|0) (1) 



(3.42) 



We now make the following Ansatz for the eigenstates of r 1 - 

\\?\...,X^ = C^\XV)C^H\?)...CV\WJ |0) . (3.43) 
These states can be related to the coefficients i^ a «--- a i in the following way : 

The state |A^\ . . . , A^) "lives" on a lattice of n sites and is thus an element of a direct product 
over n Hilbert spaces. In components it reads |A^ . . . Xn^) a , which can be directly identified 
with 

The action of r^(fi) on the states (3.43) can be evaluated with the help of the intertwiners 
(3.40) 



(3.44) 



fe=i j=i 



aw^iaw, . . . , aw) = n j n - ^i^' ■ ■ ■ . 



j" a(\f ] - m) " 



(3.45) 



fc=i i=i 



From (3.44) and (3.45) one can read off the eigenvalues of r^(fi) 
r« (Ai )|A« ...,A«> = 

_^ i JLa(u-Xi) Pr 1 A \ m m ( 3 - 46 ) 

n— ^nJr-^+n-7(T^-n^-\-) ) • 

i=i a(M - A> ; ) j-=i a( - A j W i= i a(Aj - //) i=1 y 

Inserting this expression for the special value M = A fe into (3.37), we obtain the first set of Bethe 
equations 

(a(X k )) L = Y[a(X k -\ ( t ) ) , k = l,...,n. (3.47) 
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The unwanted terms and are computed in Appendix A and their cancellation (which 
ensures that the states (3.43) are eigenstates of the transfer matrix (p.) ) leads to the following 
set of Bethe equations for the nesting 

II a ^-^ 1) )=n Ji) - P=l,.-.,n 1 . (3.48) 

i=i j=i a ( A p - A j J 

Due to our choice of grading n and n 1 can be identified as the total number of electrons and the 
number of spin down electrons respectively, i.e. n = N e = + N^, and n 1 — N^. If we define shifted 
spectral parameters according to X k = X k + ^ , we obtain the Bethe equations in their "generic" form 

feziV =n Afc 't!' § , ^=1 jv e 

^ A A (1) 4 Wi A (1) -A (1) -i 
The eigenvalues of the transfer matrix (3.35) are given by 

i=i a(fi -\\>) - =1 a l A 3 W i=1 (AJ - m) ,-=i 



(3.49) 



(3.50) 



Using the trace identities (3.19) it is possible to obtain the energy eigenvalues from the eigenvalues 
of the transfer matrix and we find 

E susy^T,T^-T-L = -2^ C0S (fc.) +2N e -L , (3.51) 

3=1 X 3 + 4 j=l 

where we have reparametrised Aj = ^coi(-j). The periodic boundary conditions (3.49) and the 
energy (3.51) are in perfect agreement with the expressions derived by Lai' 4 ' and by Schlottmann^ . 

3.5. Algebraic Bethe Ansatz with a fermionic background I (BFF grading) 
Sutherland solution 

In this section we consider a grading such that e 2 and e 3 are fermionic (representing spin 
down/spin up electrons respectively) and e 1 is bosonic (empty site). In terms of the Grassmann 
parities this means that e 2 = e 3 = 1 and e 1 = 0. We pick the reference state in the k th quantum space 
\0) k and the vacuum |0) of the complete lattice of L sites to be fermionic with all spins up, i.e. 



|0>„= , |0> = ®£ =1 |0> n (3.52) 
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This choice of grading implies that R is of the form 



R{\) = 



The L-operator is 



Z 1 














6(A) 





a(A) 








6(A) 








a(A) 





6(A) 


Q 





o 


o 




















a(A) 

















Vo 













fa 


(A) + 


&(A)e 








6(A)-a(A) 












6(A) 






a(A) 




6(A) 



-a(A) 










-a(A) 


6(A) 




(3.53) 



6(A)-a(A)/ 



b(\)el? 



b{X)el 1 
a{X)-b(X)ef 
-b(X)ef 



b{X)ef 
-b{X)ef 
a{X)-b{X)ef 



(3.54) 



The action of L k (X) on the reference state on site k is 

/ a(X) 



L k (X)\0) k = 







a(X) 



|0)» 



(3.55) 



\b(X)e k 3 -b{X)ef a(X)-b(X)J 



We partition the monodromy matrix as before in (3.24), which implies that the transfer matrix is 
now given by 

r(/x) = A n {fj.) - A^in) - D{fj) . (3.56) 



The action of the monodromy matrix on the reference state follows from (3.55) 



/(a(X)Y 



T L (X)\0) 







10) 



(3.57) 





(a(X)) L 
\ d(A) C 2 (X) (a(X)-b(X)) L J 

and inspection of (3.57) reveals that C 1 (X) and C 2 (A) are creation operators of odd and even Grass- 
mann parity respectively. We make the following Ansatz for the eigenstates of r(fi) 



\X U ..., X n \F) = C ai (X,) C a2 (A 2 ) . . . C„ n (AJ |0) 
The intertwining relations are found to be 

r(n — X) dc 

AM C C (X) = (_i)Vp+«.+<* ^—M Cp (X) am 



(3.58) 



^(-1)""- > + ^^ C ^' )A -rW 



(3.59) 



D{H) C C (X) 



1 



a(/j, — X) 



C C (X) D(fi) + — ^-C c {n)D{X) 



a(X — jj) 



C*SK) C a2 (X 2 )=r FB (X 2 -X 1 ) a a 2 1 l 1 2 C b2 (X 2 ) C^X,) 



where 



r(/i)2S =b(n)I% + a(n)n BF f d , 



=b^)i^ + a(^)n 



(3.60) 
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and H BF and U FB are the permutation matrices for the gradings e 1 — , e 2 — l and e 1 = 1 , e 2 = 
respectively. Using (3.59) we find that the diagonal elements of the monodromy matrix act on the 
states (3.58) as follows 



|A 1 ,...,A n |F> 



71 b b n 

ew """^(^n^.cA^io)^-^ 
fc =i /ai - a - f=i 



(3.61) 



(A 11 ( M )-A 22 ( M ))|A 1 ,...,AJF) = 

n 1 n 

= ( a M) L U-7— ttLKMI ) ^ (1) (M)aV.::a: f^ 

3=i a W J= i (3.62) 

+E(A t ):;:iQ t M n ^(a^io)^-^ , 

fe=l j=l 

where 

r« ( M &v/> n = (-i) £ « 4 1} (/* - aj^lI 1 !^ - K-i)l:z\lZ\ ■ ■ ■ ^ ^ ~ x 

x (-i) e =E: = " 1 l K+ i )+Er = " 1 1 %K+ i ) 



(3.63) 



Here all indices c, and c are summed over. The expression for r^(fx) is significantly different from 
the one in the FFB case treated in section 3.4., but r^(fx) can nonetheless be interpreted as the 
transfer matrix of an inhomogcneous spin model on a lattice of n sites. Our reference state |0) is 
now of fermionic nature and we have to define a new graded tensor product reflecting this fact 

(F^G) a c b d = F ab G cd (-l)^+^+^ . (3.64) 

Effectively the new graded tensor product switches even and odd Grassmann parities, i.e. e a — ► 
e a + 1. In terms of this tensor product the transfer matrix T^\fi) given by (3.63) can be obtained 
as 

r^intrX = S ^T( 1 )( M ))= S M4 1) (m-AJ^i(/"-A„-i)®...®4 1) (m-A 1 )) , 

(a(\) + b(\)ef 6(A) ef \ (3.65) 

V 6(A) ef a{\)-b{\)ef) 

In the second line of (3.65) we have explicitly written the tensor product (§> between the quantum 
spaces over the sites of the inhomogeneous model (the L-operators are of course again multiplied as 
matrices). As before F a ™ - ai must be an eigenvector of (fi), if |A X . . . A n |F) is to be an eigenstate 
of t(h). The unwanted terms can be computed in a similar way to the ones described for the FFB 
grading in Appendix B. The condition of the cancellation of the unwanted terms 

((A k )a\:X pa "- ai = ( 3 - 66 ) 
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leads to the conditions 

F^- a ^(a(-\ k )) L {rW(X k )F) an - ai , k = l,...,n . (3.67) 

To solve the nesting we first have to note that because of our change of tensor product, the L- 
operators L^(A) are not intertwined by the R-matrix r(^) defined in (3.60), but by the R-matrix 

?W& = %)<U S cd + a(n)S ad S bc (_l)«.+«c+Ve . (3.68) 

The intertwining relation 

f (A - M ) TW(\)®TW(n) = TW(»)®TW(\) f(X M ) (3.69) 

together with the choice of vacuum 

l°)i 1)= (l) - |0> C1) = ^^ijO)^ (3.70) 

can be analysed along the same lines as for the nesting in section 3.4 . It can be shown that they 
represent a model of the permutation type with BF grading (describing one species of bosons and 
one species of fermions). The resulting Bethe equations are 

m=l v 1 m ' j = l 
n 

l = Y[a(X 3 -\{ 1) ) k = l,..., ni . 

J=l 

The choice of grading implies that n and n 1 are the number of holes plus the number of spins down 
and the number of holes respectively. If we shift the spectral parameters according to 

h = \-{ > =Af>-t , (3.72) 

we obtain Sutherland's form of the periodic boundary conditions^ 4 
( A 



JJLl] = TT A J- A ^ L FT \ 2 l = l,...,N h + N i 









n 

m=l 


A, 


-A ro -i , J 








n 


*i 


~ A fc ~ 2 




_ Afe^ + 1 



(3.73) 



1= 11 J " ^ * = 1,...,^ 

j=l Aj- A fe 

The eigenvalues of the transfer matrix are 

'-/'•A, \J un/ni' If —L-,,':,,: J] 1 f — " 

J 



o(/z- A ) 



(3.74) 



4 For an odd number of lattice sites Sutherland's equation (62a) should be corrected by a factor of — 1. 
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This results in energy eigenvalues 

N H + N l , N h + N l 



L E ,-,2 1 =L-2(N h + N l )-2 J2 co < k ,) > ( 3 - 75 ) 



E 

susy 



where we have reparametrised = ^tan(-g-). 

4. A new solution of the t-J model (FBF grading) 

The third and last possibility of choosing the grading is e 1 = e 3 = 1, e 2 = 0; with e 1 representing 
spin down and e 3 spin up. This case can be analysed in precisely the same way as the BFF case so 
that we simply give the final results for the Bethe equations and eigenvalues of the transfer matrix. 
The BAE are 

= n a ( x ? i=i,...,N h +N l 

i=i 

(4.1) 

1= II o&P-lj) k = l,...,N 
j= 

We again shift the spectral parameters 



1 



\ = \-~ 2 > (4-2) 
to obtain the final expression for the new set of periodic boundary conditions for the t-J model 



L N i A (i) 

.1 2 / j=l X j - X l T 2 



(4.3) 
x (i) _ l _ i 

3=1 A k - A ] + 2 

The equivalence of the BAE (4.3) to the BAE (3.73) is demonstrated in Appendix C. The eigenvalues 
of the transfer matrix are 



Nh+Nl if / 1 \ L \ 

KMA,},^^ n ^y^ (1) ^)- (^y) j 

- (i) ^)=n-7(T^- i- n ^-^ • 

i=i a(K ~ H) \ 3=i J 



The energy eigenvalues can be obtained from (4.4) in the usual way by taking logarithmic derivatives 

n k+ n i 1 

Esusy=L- ]T ■ (4.5) 

j=l (AjJ + 4 

The string solutions of the BAE (4.3) are of a very particular structure : 

n spectral parameters X n j (j = 1 . . . n) combine with n — 1 spectral parameters 1 k (k = 1 . . . n — 1) 

16 



of the nesting into one complex string solution 

X nJ = X n + -(n+l-2j) , j = l...n 

(4.6) 

X { n \ k = X n + ^(n-2k) , fc = l...n-l . 

Note that due to the symmetry of the hamiltonian under interchange of spin-up and spin-down 
there are only three different NAB A solutions, as the other three solutions can be obtained via the 
substitution <-> Nf . 

5. Higher Conservation Laws 

In this section we derive explicit expressions for the conservation laws and (which 
involve interactions between 3 and 4 neighbouring sites respectively), using a generalisation of 
Tetel'man's method^ 22 ' 23 ' to the supersymmetric case. 

Let us define the "boost" -operator 

B = J2nH^ l+1 , (5.1) 

n 

where H^" +1 is the density of the hamiltonian given by the right hand side of (3.18). This operator 
obviously violates periodicity on the finite chain, but one can use it nonetheless in commutators 
which "differentiate" the linear n-dependence and yield formally periodic expressions. 

The integrals of motion can now be sucessively obtained by commutation with the boost- 
operator 

H {k+1) =i[B,H (k) ] 

= i[B,H {k) ] , 



where 



B = J>iJX +1 . (5.3) 



This can be seen as follows: 
If we introduce the matrix 

R(X) = ILR(A) = 6(A)II + a(A)J , (5.4) 

we can rewrite the intertwining relation (3.14) as 

R(X-ri(L n (X)®LM) = (I®LM)(L n (X)®I))R(X-») . (5.5) 

Here the L-operators are multiplied as quantum operators on both sides of (5.5). If we interchange 
the roles of matrix and quatum spaces in (5.5), we obtain the "90 degree rotated" intertwining 
relation 

Rn,n+l( X -V) { L n,l( X ) ® £„+l,l(M)) = 

(5.6) 

(4®£«+i,i(m)) (L nA (X) ® I n+1 ) R nn+1 (X-n) , 
17 



where the index "1" indicates the matrix space (which is the same for both L's, they are now 
multiplied as matrices) and n and n + 1 label the quantum spaces. The tensor product is now 
between the quantum spaces and R n n+1 is a quantum operator acting in both quantum spaces. 
From now on we will drop the matrix space index on the L-operators. In components (5.6) reads 



^((A + *)i?„,„ +1 (A))=/"'" +1 



Now we note that (5.4) implies that 

d 

! (A - (A) | - /"■" ' ' 

(5.8) 

^„ +1 (o) = ir<" +1 . 

Multiplying (5.7) by A — \i + i, differentiating with respect to A, setting A = ^, and finally acting on 
both sides of the equation with (U n ' n+1 )^ 1 '^ 1 (— l) e <*i e <*2 +c "i c ''2 from the right, we find 

[n"'"+\ L„ +1 ( M ) ® L»] = -iL n+1 (/x) ® L n (ji) + iL n+ M ® L» . (5.9) 

The tensor product is once again between quantum space and the dot denotes differentiation with 
respect to /i. From (3.18) it now follows that 

H w +1 > L n+M ® lm] = * (X+iO«) ® - £„+i(/i) ® 4.0*)) • (5.io) 

Using (5.9) it is easy to show that (up to the usual "problems" with periodicity) 

[B,r( l i)} = -it(fi) (5.11) 

and thus 

[B, ln(r (^(t(O))- 1 )] = -i^-ln (t(jm) (r(O))- 1 ) - H (2) . (5.12) 

Expanding both sides of (5.12) in powers of /x we obtain (5.2). We will now use (5.2) to obtain explicit 
expressions for higher conservation laws. According to (2.4) we can write the t-J Hamiltonian in 
terms of u(l|2) generators as 

L L L 

H susy = Hju^y 1 = - n j '' j ' +i = — k°p j ja . (5.13) 

j=i j =1 i =1 

i?( 3 ) can be obtained by commutation with the boost operator B 

H {3) =i [B,H {2) ]=i [B,H, 



= ?' \^\ffk+l,k+2 Tj-fe.fe+ll 
/ ^ L susy ' susy J 

fe=l 
L 

-i ]TiT^ K^ s f^J k _ ha J k>e J k+hS 
fe=i 



(5.14) 



Using the expressions of the w(l|2) generators (2.1) and the form of K a P — (K a/3 ) 1 given in Ap- 
pendix A it is possible to rewrite H, 3 *. in terms of fermionic creation and annihilation operators 
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L L 
j=l fe=l 

+ S't_ 1 Q] a Qj+i,-i + ^1— 1 

— 25^ S) S j+1 + (1 - Qj .! Q]+i _i + (1 - n^) q , +m 

— — — n j,-i) (5 15) 

+ Qj,-i (l-n J -_ 1 ) + Qt_ 1 _ 1 o„, 

— /i.e. . 

The u(l|2) generators Qj t<7 are given by (2.1). The next highest conservation law can be 

computed along similar lines and we find 

L 

= -2^2 K^ V K a(i K 1& f fjj e !s^Jk-l,a Jk,e Jk+l,w ^k+2,v 



k=l 

L L 

;,fe+l 

fe=l fe=l 



(5.16) 

y^pfc-i,fc+i _ 2 y^ n fc '' 
fe=l fe=l 

where p*- 1 .^ 1 i s a graded permuation operator between the sites k — 1 and fc + 1 with definition 

pk— l,fc+l j-j-fc— jjfc— 

Inspection of (5.13) reveals the expressibility of the last two terms in (5.16) in terms of u(l|2) 
generators J k a , but the physical nature of the interctions is less obvious in the resulting expression. 
The results for and given in (5.14) and (5.16) generalise trivially to the case of a lattice 
gas with an u(m\n) symmetry I 24 l . 
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A. Appendix : The Fundamental Representation of u(l|2) 

In this appendix we write down the fundamental matrix representation of w(l|2) and use it to 
compute the inverse Killing form K a/3 . We consider the fundamental matrix representation on site k, 
where we have chosen the fermionic states to be ( 1 ) T (spin down electrons) and (0 1 ) T 
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(spin up electrons) , and the bosonic state to be ( 1 ) (empty site) . The generators are given 

by 



•4,i — s k — 



J, 



fc,4 



4,7 — Qk,-1 



/o 





o\ 
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f) 


o 









/o 





u 
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(° 
















\o 





0/ 
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'k,8 



/o 
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7, — Rf — 
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V" 





1/ 



(A.l) 



The generators J k l , . . . , J k s are seen to be super-traceless, i.e. str(J) = J 33 — J 22 — J n = 0. The 
grading of the algebra is described in terms of Grassmann parities e a , which are given by 



e l = e l= e 2 = e 3 = e 8 = e 9 = 
£ 4 = e 5 = e 6 = e 7 = 1 



(A.2) 



if Q/3 = (if"' 3 ) 1 is now given by the following expression 



[14] 



K afj = str(J a J p ) 
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(A.3) 
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B. Appendix : Computation of the "unwanted terms" for the FFB grading 



In this section we compute the so-called "unwanted terms" in the expressions (3.32), (3.33), 
(3.44) and (3.45). The "unwanted terms" are characterised by containing a creation operator C with 
spectral parameter (SP) \i in place of a creation operator with SP X k (or A^ for the nesting). The 
condition of cancellation of these unwanted terms leads to the Bethe equations. In order to obtain 
the expression for A fe , we first move the creation operator with SP X k to the first place in (3.28), 
using (3.29) 

n k — 1 n 

n c ai (\) = c K (A fc ) n c K (aj n c a . s^-.-.-** , 

i=i i=i j=k+i (B.l) 

S(X k )a\...a k = r(X k _ 1 - X k ) c k k _ 1 1 a k k _ 1 »"(A fe _ 2 - ^fe)cfc_2 a fc-2 ■ • • r (Ai - A,,)^ 

To get an unwanted term, we now have to use the second term in (3.29) to move D past C bk (X k ), 
and then always the first term in (3.29) to move D (which now carries SP X k ) to the very right, 
until it hits the vacuum, on which it acts according to (3.27). This way we obtain 

Kh-») \ TT L 

The computation of A fc is more complicated. We first derive an expression for the contribution of 
A 11 (/u), which we denote by A kl . Proceeding along the same lines as in the computation of A k , we 
find 



L = i «( A fe - A;) 

x KAfc-A^; r(A fc -A 2 )££ . . . r(X k - X k _,) x (B.3) 

x r(X k - Afc+i)^^ 1 ^ r(X k - K+2) bl + + ld k +2 • ■ • r (A fe - A„) fc ™ d ^ x 
x (a{X k )) L 6 dn iA (-I)"" 1 . 

The (5 d i j stems from the action of A 1 d i (A fc ) (which is what we get after moving A past all the 
C's) on the vacuum. We also had to include a <5 6fc x due to the fact that in (3.33) we denoted by b k 
the index of the C with SP /1 (this means that on the l.h.s. of (B.3) b k is really equal to 1, too). The 
contribution of A 22 (p) differs only in factors 5 b ^ 2 <5 d ± 2 instead of 5 bk 1 S d i 1; so that the result 
for A fc = A fc l + A k 2 is found to be 



))<- X 



(B.4) 



x KA fe -A 1 )^r(A fe -A 2 )^ ■ • • r(X k - Afc_i)&*~*d*~* (-*)" 
x r(A fe - A fc+1 ) b ^ + ^+ fc 1 _ i r(X k - A fe+2 ) fo ^+ 2 1 d ^ +2 . . . r(A fe - A„)^™_ 2 

This expression can (and must be) simplified by carrying out the contractions over the summation 
indices c l5 . . . , c k . Noting that 

KAi - A 2 )^ r(A 2 - XJUll = = S aiCi S a2C2 , (B.5) 
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we are able to perform all c r summations with the result 

fc-i 



s(\ k y a \:% Hx k ■ ■ ■ r ^ x ^KSt\ = W ■ (b.6) 



•••rlAfc - A k-l)b k _[d k _' 2 - 1 

Now we transform the remaining r-matrices into L-operators, using the identity 



r(\)f d = (r(A)nW)^ (-1)^ = -(r(\)Tl {1) ) a J = -L^(X) a J . (B.7) 

The second equality holds, because e x = e 2 = 1. Thus we obtain our final form for the unwanted 
terms due to A u ([i) + A 22 (^i) 



{AkF)b ,..K = fcM JJ {a{Xk))L F a n .., kh . 1 ., 1 



X 



<*(/*- Afc) a{\ k -\ t ) 

Ifk (B.8) 



We now insert (B.8) and (B.2) into the condition (3.36) for the cancellation of the unwanted terms, 
and multiply the resulting equation by the inverse of S(\ k )a\'.'.' b a k , which satisfies (iS _1 (A fc )) j 1 '"^* x 
( 5 '( A fe))a 1 1 '.'.'.a fe = IL=i <*a i)Pi > and which is computed via (B.5). 



After some trivial rearrangements we arrive at 

7=i a ( x i- x k) 

i=fk 



- L n H X k - X n)b n a n ^ ' ■ ' L k+l( X k ~ X k+l) b a 



Mn-2 r (l) 



> Pl b k r y- 1 ) (B.9) 

= L l-l( A fc - A '=-l)p^_ fc 1 6fc_ 1 - £ 'fe-2( A fc _ Xk -^pS 2 bS 2 ■ ■ ■ L i ^ X k - A l)p 2 1 hi X 



This implies (3.37). 

The unwanted terms (3.44) for the nesting can be computed along similar lines and we easily 
find that 

j^k 

(B.10) 

k a(AW- M )M^ 1) -AS 1) )iia(A l -A«) ' 
The cancellation of the unwanted terms for the nesting implies (3.48). 
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C. Appendix : Equivalence of the BAE 



In this appendix we establish the equivalence of our new set of BAE (4.3) with the set obtained 
by Sutherland (3.73). We use a method developed by Bares et.al. in [25] (see Appendix C of their 
paper) to show the equivalence of the solutions of Sutherland and Lai (see also [26] ) . 

We start by expressing the second set of our new BAE 



' x j- * \ L Nl A (1) X i 



N h +N l (1) _ ~ { 

n A k A 3 2 . _, N 

.(1) r ,i k-L,...,J\ l 



(C.l) 



X w — X 4- 
3=1 A k A ] + 2 



as a polynomial equation of degree N h + 

N h +N l . N h +N l 

P (w)= n {™-\-\)~ n (^-a,+|)=o . (c.2) 

j=i j=i 
Among the roots Wj, j — 1 . . . N h + of (C.2) we consider roots w 1 ,. . . , w N ^ , which we identify 
with . . . , Al^. The N h other roots of (C.2) we denote by w' •. Using the residue theorem we 
can derive the following equality 

where Cj is a small contour around Wj . The branch cut of the logarithm extends from z n = A ; + | 
to ^ p = A ; — |. By deforming the contours on the r.h.s. of (C.3) we arrive the the following equality 

where the last term on the r.h.s. comes from integration around the branch cut. The form of the 
polynomial p now implies that 

P( z n) = - II { A l- A m + i) 
m—1 

P(Zp)= LI ( A l-~ X m-i) ■ 
m—1 

Inserting (C.5) into (C.4) and exponentiating the result we obtain the identity 

"j. A ; -Wj + % = ^ A; -w' k -i Nh J^ 1 A ; -A m + z g 

Now we use (C.6) in the first set of BAE in (C.l) with the result 

hill) = TT ^rn + i T[ X J- w k - 1 2 l = 1 ,..., Nh + N . (c.7) 

m^l 
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(C.5) 



This is precisely the first set of BAE in the Sutherland form (3.73), if we make the identification 
w' k = . The second set of the BAE (3.73) is also fullfiUcd by the spectral parameters A^ 1 ', 
because they are roots of the equivalent polynomial equation (C.2) (this is because the BAE of the 
nesting for the Sutherland solution and the new solution are identical). Thus we have established 
the equivalence of the BAE (4.3) and (3.73). 
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